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ABSTRACT 

The bimodality in galaxy properties has been observed at low and high redshift, with a clear 
distinction between star-forming galaxies in the blue cloud and passively evolving objects in the red 
sequence; the absence of galaxies with intermediate properties indicates that the quenching of star 
formation and subsequent transition between populations must happen rapidly. In this paper, we 
present a study of over 100 transiting galaxies in the so-called "green valley" at intermediate redshifts 
(z ~ 0.8). By using very deep spectroscopy with the DEIMOS instrument at the Keck telescope we 
are able to infer the star formation histories of these objects and measure the stellar mass flux density 
transiting from the blue cloud to the red sequence when the universe was half its current age. Our 
results indicate that the process happened more rapidly and for more massive galaxies in the past, 
suggesting a top-down scenario in which the massive end of the red sequence is forming first. This 
represent another aspect of downsizing, with the mass flux density moving towards smaller galaxies 
in recent times. 

Subject headings: galaxies: evolution, galaxies: luminosity function 



1. INTRODUCTION 

The distinction between blue star-forming galaxies and 
red, passively evolving ones has been k now for a long 
time. In his seminal classification article, iHubbld (|1926[ ) 
classified galaxies into two main groups, spirals and ell- 
lipticals. Although unintended as a timeline for secular 
evolution - in Hubble's own words 'The entire classifica- 
tion is purely empirical and without prejudice to theories 
of evolution' (|Hubbldll927f) - ellipticals were referred to 
as 'early- type galax ies' and spirals as 'late types' (see also 
iFortson et aHl2011l) . 

Our knowledge of stellar populations in galaxies has 
evolved considerably since then. Today we know that 
spiral galaxies are blue due to young stars, and that their 
current star formation rates are higher th an in ellipticals, 
with a younger population, on average. IDavis fc Geilirl 
(|1976() were the first to show that clustering properties of 
galaxies depend on galaxy type, with ellipticals showing 
stronger clustering properties than spirals. We also know 
that elliptical galaxies in the local universe are preferen- 
tially found in environments where the galaxy number 
densi ty is higher such as clusteres of galaxies (jDresslerl 
1980), which gives us important clues regarding the late 
evolution of galaxies. The same relation has been shown 
to extend to groups of galaxies ([Postman fe Gellerlll984fl 
and e nvironments such as superclusters (fGiovanelli et al.1 
1986). This leads to the conclusion that spiral galaxies 
somehow transform into ellipticals at later times, for ex- 
ample through secular evolution or mergers with other 
galaxies. However, the existing bimodality in the galaxy 
distribution, with a clear distinction between blue spi- 



rals and red spheroids, and its establishment still remain 
a puzzle. 

Other works have correlated the bimodality in galaxy 
prope rties with galaxy colors. IVisvanathan fc Sandagd 
were the first to show the correlation between 
colors and magnitudes of cluster galaxies in the local 
universe. The bimodality in galaxy colors and their 
correlations wit h galax y morphology was first noted by 
iTakamiya et al.l ([1995). More recently, the color bi- 
modality has been observed at low reds hifts (z ~ 0.1) 
with t he Sloan Digital Sky Survey (SDSS: IStrateva et al.1 
120011: iBaldrv etaLl l2004fi . Later, many authors have 
shown this bimodality extends to ear li er times, at 
least out to z ~ 1.0 (IK00 et~all Il996t iBalogh et al.l 



19981: Urn et al. 



Willmer et al.l 



20021 : IBell et alJl2004t IWemer et a' . 2005; 
2006). As is the case for galaxy mor- 
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phologies, several recent works have shown that galaxy 
colors are correlated with environme nt density (e.g., 
iBlanton et aLll2005t iCooper et aT]|2006l ). If the color dis- 
tribution of galaxies can be described as two distinctive 
peaks, then one can denominate the minimum at inter- 
mediate color values the "green valley". The question 
then arises: why does such a minimum exist, instead of 
a homogeneous distribution across the color-magnitude 
diagram? 

In order to study this intermediate population, 
one needs t o de termine proper selection criteria. 
iWyder et al.l (|2007t ) have measured number densities in 
the color-magnitude diagram at low redshift (z ~ 0.1) 
and found that the GALEX NUV — r color represents 
more efficient criterium to select green valley galaxies. 
Ultraviolet emission originates from recent (over the last 
tens of millions of years) star-forming regions in the 
galaxy, while the r band is more sensitive to the bulk 
of stellar mass, formed over the course of the galaxy's 
history. By subtracting the two, one can clearly distin- 
guish two populations, one actively forming stars and 
the other older, more passive, with a dynamic range of 
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about six magnitudes in color. Nevertheless, we need 
to be careful when defining transition galaxies; the per- 
centage of star-forming, obscured galaxies is high. These 
are typically galaxies with high ratio of current to past 
star formation, b ut which present redder colors due to 
dust obscuration ([Martin et al.1 120071 : ISalim et al.ll2009L 
see also section I3~4l . 

The low number density (or number deficit) of galax- 
ies in the green valley indicates that the transition be- 
tween both groups occurs rapidly. A number of works 
attempt to explain why there is such a rapid evolution. 
IMenci et ail ([20051) postulate that supernova winds and 
outflows can play an important part, driving out gas 
that f uels star formation in the galaxy. iDi Matteo et al.l 
(2005), on the other hand, have produced a hydrody- 
namical simulation of a major merger event which shows 
that, after a period of brief increase of star formation, 
the supermassive black hole created in the center of 
the merger remnant drives st rong outflow winds that 
rapidly quench star form ation. iNandra et al.l (|2007T) and 
ISchawinski et al.l f|2009T ) provide observational support to 
this hypothesis, finding a large number of active galac- 
tic nuclei (AGN) in the green valley and concluding that 
feedback from such objects might be somehow quench- 
ing tfifi^tai^rmation process. This is further supported 
by iCoil et ail (|2008l) , who note that the coadded spec- 
tra of green valley galaxies at z ~ 1 are distinct from 
coadded spectra of blue and red ga laxies, with l ine ra tios 
that show increased AGN activity. IPeng et "ail ([2010D go 
further and state that quenching is a combination of a 
mass-dependent quenching mechanism, proportional to 
the galaxy's star formation rate, and the effect of the hi- 
erarchical assembly of dark matter haloes. Nevertheless, 
IMendez etail (|2011[ ) note that the merger fraction in the 
green valley is low, from the quantitative morphological 
analysis of 0.4 < z < 1.2 optically selected transition 
galaxies. The authors then conclude that mild external 
processes (such as galaxy harassment) or secular evolu- 
tion are the dominant factors at these redshifts. 

This rationale assumes, however, that galaxies move 
in one direction only in the color-magnitude diagram, 
from blue to red. Nevertheless, there are several 
cases in the literature that indicate ongoing and/or 
recent gas accretion, in addition to comparatively 
high star form ation rates in elliptical galaxies (e.g., 
Morganti et al.ll2006tlKayirai et al.ll2008HCatinella et al.1 
20101 : ISchiminovich et ail 120101) . Although these cases 
are interesting, they are relatively rare, and one can in 
general assume that galaxies with intermediate colors are 
currently in the process of quenching their star forma- 
tion. 

To infer how rapidly galaxie s are moving across the 
green valley. IMartin et ail (|2007D (hereafter Paper I) have 
used spectroscopic features in green valley galaxies to 
obtain information on their star formation histories (see 
section [5] for details) . Along with measurements of typ- 
ical galaxy masses and number densities in the color- 
magnitude diagram, the authors have been able to de- 
termine the mass flux across the green valley at low red- 
shifts (z ~ 0.1). The measured value of p — 0.033 M Q 
yr _1 Mpc~ 3 agrees remarkably well with expectations 
from the growth of the red sequence and the d epletion of 
galaxies in the blue sequence at such redshifts ■: Hell et al.l 



120041: iFaber et al.l [2007|), despite the shortcomings of the 
simple model assumed. Furthermore, Paper I presents 
evidence for an increase in number density of AGN in 
the green valley (from measured [OIII] luminosities), sup- 
porting the aforementioned studies that relate AGN ac- 
tivity with the quenching of star formation (although 
the correlation between [OIII] luminosities and quench- 
ing timescales is not unequivocal). 

In this work, we apply the methods introduced in pa- 
per I to a sample of galaxies at higher redshift. This 
will allow us to measure the mass flux in the green valley 
at intermediate redshifts, comparing with results found 
for galaxies in the low-redshift universe. Although the 
technique is simple, the data acquisition process proves 
challenging, since it requires reliable measurements of ab- 
sorption features in galaxies at redshifts of (z ~ 0.8). We 
have set out to obtain the required data, and the work 
presented here includes the deepest spectra of green val- 
ley galaxies to date. 

This paper is organized as follows. In section [2] we 
describe in detail the methodology used to infer star for- 
mation histories of galaxies, including modelling of stellar 
populations. In section [3] we describe the observations 
and data reduction used to produce the spectra which 
were then used to measure star formation histories and 
quenching timescales, in addition to ancillary data used 
to measure number densities in the color magnitude di- 
agram, luminosity functions and extinction correction. 
Section [4] shows our results, including the measured mass 
flux density at z ~ 0.8 and in section [5] we discuss those 
results in light of galaxy evolution models. We summa- 
rize our findings in section [51 Throughout this paper, 
we use AB magnitudes, and assume standard cosmology, 
with H = 70 km s 1 Mpc 1 , ti m = 0.30, and Q A = 0.70. 

2. METHODOLOGY 

The method to study the mass flux i n the green valley 
has been introduced by IMartin et al.l (|2007[ ) . We sum- 
marize here the description presented in that work. 

2.1. The mass flux density in the color magnitude 
diagram 

In order to measure the mass flux for a given color 
in the color-magnitude diagram, M{r, y), where y is the 
NUV — r color of the galaxy, we can assume that 



M(r,y) =< M(r,y) > x 



dy 
dV 



(1) 



where < M(r,y) > is the average mass of a galaxy in 
that color- magnitude bin, and dy/dt is how fast galaxies 
are moving through the same bin. In practice, we mea- 
sure average values in a two-dimensional bin. Dividing 
the above equation by the comoving density probed, we 
calculate 

p{r,y) = $(r,y)< M(r,y)>^. (2) 

In this case, the comoving mass density p(r, y) is simply 
the comoving number density ($(r, y)) multiplied by the 
typical galaxy mass in that bin (< M(r,y) >). What 
we propose to measure is the mass flux density p. Since 
we can constrain <f>(r, y) and < M{r,y) > independently 
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from the star formation histories, we are only left with 
the task of measuring the color evolution rate for a given 
galaxy or bin. 

2.2. Star formation histories 

In order to measure the color evolution rates, we make 
some simplifying assumptions: 

1. Galaxies only move towards redder colors, i.e. we 
do not consider starbursting red galaxies moving 
downward in the color-magnitude diagram; 

2. The star formation histories in all galaxies are de- 
scribed by a constant star formation rate for the 
first few Gyr (approximately 6 Gyr) followed by 
a period of exponentially declining star formation 
rates: 



In Figure [TJ we show 5 tracks on the _Dn(4000) vs. 
Hs,a diagram given by 5 distinct 7 values (0.5, 1.0, 2.0, 
5.0 and 20.0 Gyr" 1 ). These m odels were prod uced with 
the IBruzual fc Chariot! ([20031 ) models - with IChabrierl 
(2003) initial mass functions, Padova 1994 stellar evolu- 
tionary tracks and solar metallicities - by varying ages 
throughout star formation histories as described in equa- 
tion [3] The range of values for to is limited to to < t z , 
so that the oldest models are always below the universe 
age at the highest redshift we measure. This does not 
have a major impact on our results, since to is typically 
larger than the time when Z?„(4000) and H$a have sta- 
bilized. We notice very distinct tracks, with the strongly 
quenched models (higher 7 presenting higher Hs,a val- 
ues). 



SFR(f) 



SFR(t = 0), 
SFR(i = 0)e- 7t , 



t < t 

t > t 



(3) 



We discuss the implications of this model in Section 15.21 
In order to measure the exponentia l index 7, we ap- 
ply the same methodology described in lKauffmann etaLl 
(2003) to measure the rest-frame 4000 A break and the 
equivalent width of Hg absorption. The first is created by 
the accumulation of a large number of ionized metal ab- 
sorption lines shortward of 4000 A (Bruzual 1983). The 
elements responsible for these lines are multiply ionized 
in hot stars, resulting in smaller opacities in the inte- 
grated spectra of younger galaxies, and thus the break 
correlates with the age of the stellar population. The 
latter is mainly present in the stellar photosphere of 
smaller early- type stars (mostly A stars), which domi- 
nate 0.1 — 1.0 Gyr after a starburst event (and after the 
O and B stars evolved off the main sequence). 

In this work, we use the same definitions as in Paper I, 
in order to maintain uniformity across different redshifts. 
The first index, D„(4000), is defined as the ratio of the 
average flux density F v in the bands 3850-3950 and 4000 - 
4100 A, following the definition of iBalogh etaLl (|1999f ). 
The narrow bands ensure that the ratio is weakly depen- 
dent on flux calibrations and other broad band effects 
that may arise from data reduction. For the latter, H$,a 
is the absorption equivalent width; the continuum is de- 
fined by fitting a straight line through the average flux 
density between 4041.60 and 4079.75 A, on one end, and 
4128.50 and 4161.00 A, on the other. The equivalent 
width, then, is given simply by 



H, 



8, A 



E 



A=4083.5 



F„ 



(4) 



Kauffma nn et al.l f|2003T ) showed that these indices 
trace a well defined region in a Dn(4000) vs. H$,a di- 
agram. Furthermore, different star formation history 
tracks - for instance a single starburst event versus a con- 
tinuous star formation rate - trace distinct regions within 
this diagram. We therefore use our measured spectral 
indices as defined above to distinguish between different 
star formation histories, in our case as given by different 
7 values. 



< 

< 

to 

PC 




1.0 1.2 1.4 1.6 1.8 
Dn(4000) 
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Fig. 1. — Evolution of D n (4000) vs. Hg^ for models with differ- 
ent star formation histories. Solid, dotted, dashed, dotted-dashed 
and triple-dotted-dashed lines (from left to right) represent 7's 
equal to 0.5, 1.0, 2.0, 5.0 and 20.0 Gyr -1 , respectively. Therefore, 
models to the right are quenching star formation more rapidly. 

Finally, the measurement of dy/dt is given by the width 
of the color range in magnitudes used to determine the 
green valley divided by the duration of time when the 
galaxy is within those boundaries. We illustrate this in 
Figure [21 where we plot the NUV — r color as a function 
of time after the initial period of constant star formation, 
with the green valley limits indicated by horizontal lines. 
As expected, models with higher values of gamma (i.e. 
that are quenching their star formation more rapidly) 
change colors faster. The timescales to cross the green 
valley in these models varies between 0.26 and 2.7 Gyr. 

3. SAMPLE, OBSERVATIONS AND DATA PROCESSING 

3.1. Sample selection 

We have selected a sample of 163 green valley galaxies 
from surveys at intermediate redshift surveys. Preference 
has been given to galaxies in the Extended Groth Strip 
(EGS), which is the target of the ongoing All- wavelength 
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Fig. 2. — NUV — r colors of model galaxies as a function of 
time after the onset of exponential decline in SFR. Line types are 
the same as in Figure [l] The thick horizontal lines indicate the 
threshold we have used to determine the green valley, dy/dt will 
then simply be the width of the magnitude band in colors divided 
by the time take for a galaxy to cross both thresholds. 

Extended Gr oth Strip Internat ional Survey collabora- 
tion ( AEGIS; iDavis et al.l [20071) . This field is optimum 
for galaxy evolution studies, since the AEGIS collabora- 
tion has produced large amount of ancillary multiwave- 
length data. In particular, a large subsample of galax- 
ies has spectro scopically determined redshfits w ith the 
DEEP2 survey pavisll2003t iNewman et~aLlll968[ ). How- 
ever, in order to facilitate observations throughout the 
year, we have also used data from the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLSjB, from all 
deep fields observed. 

All galaxies observed in this work have CFHTLS pho- 
tometry in all 5 bands (u,g,r,i,z). This is required 
to properly constrain SDSS r' band magnitudes in the 
galaxy's restframe, while simultaneously yielding accu- 
rate GALEX NUV restframe magnitudes. In order to 
match the magnitude limit of the spectroscopic sample 
in DEEP2, we have determined a threshold of M r < 24 
for all objects in our sample. The CFHTLS Deep Field 
D3 overlaps with the EGS, and in that case we match 
CFHTLS and DEEP2 sources. We preselect galaxies 
with redshifts between 0.55 < z < 0.9; DEEP2 rcdshifts 
were used for the selection in the EGS field, CFHTLS 
photometric redshifts are used otherwise (approximately 
50% of our final sample). 

The final sample of 163 observed galaxies was selected 
at random from the available targets, in an attempt to 
maximize object count in a given DEIMOS mask. Not 
all objects in a given area were available for observation 
in a single pointing, given the number density of green 
valley galaxies and the design constraints for slits in each 
mask. Preference was given to objects in the green valley 
after extinction correction, but no additional criteria are 

4 http://www.cfht.hawaii.edu/Science/CFHTLS/ 



used in source selection. 

All sources have been k-corrected to redshift z = 
using the Kcorrect code (version 4_2; lBlanton fc Roweid 
2007), so that we could select them based on rest-frame 
NUV-r colors. In Figure [3] we show the result of the k- 
correction processes, showing all galaxies below the com- 
pleteness limits in CFHTLS-Deep. We also indicate the 
color selection criterium, showing the color limits for the 
green valley galaxies as dashed lines. 




Fig. 3. — Color-magnitude diagram of all CFHTLS sources be- 
tween 0.55 < z < 0.9 in logarithmic contours. Two distinct popula- 
tions are easily distinguishable, with the "green valley" in-between. 
The dashed lines indicate the color selection criterium used in this 
work. 



3.2. Observations and Data Reduction 

Our observations were undertaken with the DEIMO S 
instrument on the Keck II telescope (Faber 2003). 
DEIMOS is a multi-object optical spectrograph, which 
allows up to nearly 100 objects per mask in a field of view 
of 16.7 arcmin x 5.0 arcmin. The number of GV galax- 
ies per mask oscillated between 20-30 obj/mask; other 
galaxies in the aforementioned surveys were then used 
to complement observations. We used the 1200 mm -1 
grating for a resolution of R ~ 5000. A grating angle 
centered on 7500 A was chosen in order to cover all the 
necessary wavelength range needed to measure D„ (4000) 
and Us. a at the redshift range of our sample. In prepar- 
ing the masks, preference was given to objects in the 
extinction-corrected green valley (section I3.4[) . 

In order to achieve the required signal-to-noise level 
(S/N ~ 3 per pixel, or S/N ~ 2 for galaxies fainter than 
r > 23) , we exposed for a total of up to 9 hours per mask, 
weather and sky availability permitting. This comprises 
the deepest spectra of intermediate redshift green valley 
galaxies taken to date. The observations are summarized 
in Table [1] 

Preliminary data reduction was performed with the 
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DEEP2 pipeline^- The pipeline is currently optimized for 
measuring redshifts in galaxies observed with the DEEP2 
survey; therefore, we have used it to select individual 
spectra, rectify slits and subtract the sky background. 
Extraction to one-dimensional spectra, refinement of the 
wavelength solution through the use of sky lines and fur- 
ther data analysis have been done with custom IDL pro- 
cedures. All redshifts have been remeasured to ensure 
correct wavelengths for measuring D„(4000) and H^. 

When measuring D„(4000) and Us, a values for indi- 
vidual galaxies, we have resampled the spectra so that 
the wavelength res o lution per pixel is the same as the 
iBruzual fc Chariot! (pOOl models, i.e. 1 A. In all cases 
we mask out regions contaminated by sky lines. We show 
the resulting spectra in Figure [4] 




3600 3800 4000 4200 



TABLE 1 
Summary of DEIMOS Observations 



Mask 


Field 


Observing 


Total integration 


Number of 






date (UT) 


time 


GV galaxies 


AG0801 


EGS 


May 2008 


32400 


45 


AG0803 


EGS 


May 2008 


25200 


48 


DZLE01 


DEEP2 Field 2 


May 2008 


3600 


18 


D401 


CFHTLS D4 


August 2008 


16200 


37 


XMM02 


XMM-LSS 


August 2008 


7200 


12 


AG0901 


EGS 


April 2009 


28800 


33 


AG0902 


EGS 


April 2009 


28800 


30 



in distinct regions of the sky (Dl, D2, D3, D4), the in- 
fluence of cosmic variance is minimized. 

In determinin g number den sities, we have used the 
l/Vmax method (| Schmidt! 1 19681 ). This takes into account 
the magnitude limits of the survey and the potential of 
missing low-luminosity galaxies. In that sense, we per- 
form a k-correction of every source in the survey between 
the limiting redshifts (0.55 < z < 0.9) to determine the 
maximum distance at which we would be able to detect 
it (below the z = 0.9 cutoff), taking into account all 5 
bands used to select sources. We denote the maximum 
redshift as 




3600 3800 4000 4200 

X (A) 



Fig. 4. — Spectra of three green valley galaxies plotted in the 
region of interest (between 3600A < A < 4200A). Observed r 
magnitudes are approximately 21.5, 22.5 and 23.5, respectively. 
Spectra have been re-sampled to match the model resolution of 
1 A per pixel. We also highlight the position of [Oil], Call H & 
K and lines. In all cases the calcium and lines are clearly 
distinguishable. 



3.3. Number densities and the luminosity funcion 

In order to determine the number density per bin in the 
color-magnitude diagram, we have used only CFHTLS 
photometric data, including photometric redshifts when 
spectroscopic data are not available. These redshifts are 
precise enough for this exercise, and the increase in sam- 
ple size and area greatly improves the statistics. Fur- 
thermore, since the CFHTLS includes 4 different fields 

5 The pipeline is available at http://deep.berkeley.edu/spec2d/ 



0.9). (5) 

Zmin is simply 0.55, since there is no brightest magnitude 
cutoff. The maximum volume for each galaxy will be 

Kmax - 3 Vl80J U + W) 3 (l + z mi „) 3 J' (j 

where A is the angular area in the sky in square degrees 
occupied by the four fields of the survey and Dl{z) is 
the luminosity distance to redsfhit z. 

Finally, the number density in the color-magnitude di- 
agram is 

in units of Mpc~ 3 mag -2 , where AM r and Ay are the 
magnitude and color bin size, respectively. We show our 
results in Figure [3] The colored regions indicate number 
densities as measured from the CFHTLS sample, while 
the black dots are the objects targeted in this work. 
Two regions of higher number density stand out, evi- 
dence of the color bimodality still present at these highe r 
redshifts, similar to the results of iWillmer et al.l (|2006f) . 
Howeve r, a comparison w ith similar work at redshift 
z ~ 0.1 (jWvder et a l. 2007) shows that the number den- 
sities in the red sequence are smaller, showing there has 
been a significant growth since z ~ 1 dFaber et al.l [20071 
see al so Figures [3] and [5] and Figure 7 in iWvder et all 
(2007)). We discuss this in more detail in section [5] 

We generate luminosity functions by multiplying the 
number densities by the color bin size Ay. In Figure [51 
we show the resulting luminosity functions for three dif- 
ferent color bins (NUV - r = 1.75, 3.75, 4.75), and com- 
pare those w i th th e same functions at z — 0.1, from 
IWvder et al.l (|2007l) . Although galaxies have been k- 
corrected to different redshifts (z = for our sample 
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Fig. 5. — Diagram of number density of galaxies in each two- 
dimensional bin in the color-magnitude diagram, in #/Mpc — 3 . 
Two populations are distinguishable: the red sequence at NUV — 
r ~ 5 and the blue sequence at NUV — r ~ 2. Green valley galaxies 
used in this study (black points) fall in-between these two popula- 
tions, at 3.5 < NUV - r < 4.5. 



$(M) = 0.41n(10)$*10- a4 ( M - M *» Q+1 )exp[-10- a4 ( M - M *)] 
where M* and $* are the free parameters in the fit. 




and z = 0.1 for lWvder et all (|2007D ). the typical shift in 
magnitude is M r =0.1, and thus this difference has no 
impact on our results. 

Assuming the dynamic range is insufficient to con- 
strain the faint-end slope of the luminosity funcion, we 
fix it at the same values as for the lum inosity fun c tion a t 
low redshift, followin g th e strategy of [Bell et al.l (|2004l) . 
iWilhner et ail (|2006l ) and lFaber et al.l (|2007l ). As Figure 
[6] shows, this is satisfactory for the blue sequence and 
the green valley; the luminosity function of the red se- 
quence, however, clearly shows a deficit of faint galaxies 
when compared to its low-redshift counterpart. In this 
case we also fit a Schechter luminosity function with the 
faint-end slope as a free parameter. We discuss this evo- 
lution in section [5] 

We see that (1) green valley and red sequence num- 
ber densities are smaller than in the local universe, while 
number densities are similar for the blue sequence, and 
(2) all luminosity functions are shifted towards higher 
luminosities - a fact that has been extensively observed 
at high redshift and has been cited numero usly as ev- 
idence for another mode of downsizing (e .g iBell et al~l 
l200i IBundv et al.ll2f)M iFaber et al.ll20t)l . It should 
be noted that the red peak in the color-magnitude oc- 
curs at slightly different colors at low and high redsfhit, 
with the red sequence being slightly bluer at z ~ 0.8. 
However, the variation between the NUV — r = 4.75 
and NUV — r = 5.75 luminosity functions at z = 0.1 
are small, and chosing a different color range would not 
change our results. In Table [5] we present the resulting 
parameters from a Schechter function fit in each case, as 
well a s the low-redshift values presented in lWvder et al.l 
(|2007[ h The luminosity function is described as follows: 



Fig. 6. — Luminosity functions at 2 ~ 0.8 measured from 
the CFHTLS Deep survey. Blue represents the blue sequence 
(NUV - r = 1.75), red, the red sequence(NUV - r = 4.75), and 
green, the green valley (NUV — r = 3.75). Solid lines represent the 
fit for a fixed faint-end slope equal to values found at low redshift. 
Additionally, the red dotted line represents the Schechter function 
fit with the faint-end slope as a free parameter. The dash-dotted 
dark-green line and symbols represent the extinction-corrected lu- 
minosity function for green- valley galaxies. Error bars are given 
by Poisson statistics only; cosmic variance errors should be mini- 
mal, since we use all four deep fields in CFHTLS to compute the 
luminosity functions. The das hed lines indicate luminosity func- 
tions measured at z = 0.1 bv lWvder et al.l 1120071) . for comparison. 
All functions at higher redshift are shifted towards brighter mag- 
nitudes. The blue sequence shows similar number densities, while 
red and green galaxies are rarer than in the local universe. 



3.4. Extinction correction 

Extinction correction is of fundamental importance to 
identify quenching galaxies in the green valley. This is 
mainly due to contamination by heavily obscured star- 
forming galaxies. This is already an important problem 
at z ~ 0, but we expect it to be worse at higher red- 
shifts, since the number density of LIRGs and ULIRGs is 
found to increase towards earlier times ([Le Floch et al.1 
120051: iMagnelli et "all 120091 ). To illustrate the issue, we 
show in Figure [7] the fraction of 24 /mi detected galaxies 
in the Groth Strip per two-dimensional bin in the color 
magnitude diagra m (for a description of the dataset, see 
IDavis et al.1 12007T ) . We can see that most of the green 
valley galaxies on the bright end (up to ~ 65% in some 
color-magnitude bins) are also 24 /im sources, indicative 
of dusty star forming galaxies, instead of the quenching 
objects for which we are searching. 

In an attempt to decrease contamination from such 
galaxies, we used an independent SED fitting result 



Mass flux in the green valley 7 



TABLE 2 

SCHECHTER FUNCTION PARAMETERS 



Sample 




Log 4>* 


M» 


a 


Blue sequence (NUV — r 


= 1.75) 


-3.05 ±0.05 


-21.73 ±0.05 


-1.465 


Red sequence (NUV — r 


= 4.75) 


-3.39 ±0.05 


-22.90 ±0.06 


-0.579 


Red sequence (free a) 




-3.17 ±0.05 


-22.10 ±0.11 


0.24 ±0.14 


Green valley (NUV — r - 


= 3.75) 


-3.26 ±0.05 


-22.20 ±0.06 


-0.357 


Green valley (extinction- 


corrected) 


-3.60 ±0.07 


-23.07 ±0.08 


-0.357 


Blue sequence (z = 0.1) 




-2.871 


-20.331 


-1.465 


Red sequence (z = 0.1) 




-2.962 


-20.874 


-0.579 


Green valley (z = 0.1) 




-2.775 


-20.711 


-0.357 
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Fig. 7. — Fraction of 24 fim detections as a function of color and 
magnitude. The thick horizontal lines indication the region of the 
color-magnitude diagram defined as the green valley in this work. 
Up to 65% of the green valley galaxies in a given bin are probably 
dusty starbursts, and not galaxies in the process of quenching star 
formation. 

for a number of g alaxies in the Extended Groth Strip 
(|Salim et a l. 2009), where extinction is a free parameter. 
This sample of ~6000 objects was then used to calculate 
the number densities in the color-magnitude diagram the 
same way as described in section 13.31 The result can be 
found in Figure [SJ The main difference when compar- 
ing this diagram with Figure [5] is a decrease in num- 
ber density in the green valley and a clearer distinction 
between the blue and red sequences, as expected. We 
have also plotted the corresponding luminosity function 
for extinction-corrected green valley galaxies in Figure [SJ 
There is a decrease in $*, as expected; in addition, we 
notice a brighter value of M*, since this represents the 
magnitude of galaxies corrected for internal extinction. 

4. RESULTS 

We have measured D„(4000) and Hs,a indices for all 
green valley galaxies observed. A number of galaxies pre- 
sented unrealistic values - outside the expected range of 
1.0 < L>„(4000) < 2.0 and -4.0 < H S)A < 8.0. These 
objects present a number of issues, including signal-to- 
noise ratios S/N < 1 per pixel, or they are too close to 
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Fig. 8. — Same as Figure[5] after app lying extinction correction 
of absolute magnitudes (see section [3.41 1. 

the edge of a mask and a proper extraction was not pos- 
sible. Excluding these, we are left with 105 green valley 
galaxies. These objects are shown in Figure|Hl along with 
the models presented in Figure [TJ The observed values 
agree well with the models, following the trend of lower 
Hg,A values for higher L>„(4000). 

In selecting galaxies for the extinction-corrected green 
valley, our sample was reduced, due to smaller num- 
ber densities and only a fraction (approximately 50%) 
of galaxies in DEEP2 (and no galaxies in the remaining 
CFHTLS fields) having reliable SED extinction measure- 
ments. For this reason, we have expanded the color range 
to include 3.0 < (NUV — r) cxt < 4.5 colors, in an attempt 
to improve our statistics. Extinction-corrected galaxies 
are highlighted in Figure [3] as red symbols. 

Figure [5] also shows the D„(4000) and values for 
models of green valley galaxies defined at NUV — r = 4.0, 
the median color value of our green valley sample. The 
straight diagonal lines indicate equidistant lines from two 
adjacent models, in which distance is defined by normal- 
izing the Hs,a index as Us, An — Hs,a/12 to reflect the 
expected dynamic range in that measurement. 

In Figure [10] we compare the quenching timescales we 
find with the obtained values at z = 0.1 in Paper I, 
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Fig. 9. — D n (4000) and a indices of green valley galaxies with 
errorbars. The different lines on the plot indicate the expected 
values from models with varying 7 values, as in Figure [T] Black 
points are all galaxies that fall in the green valley region, prior to 
making a correction for extinction in the selection. Red symbols 
indicate galaxies in the green valley after the extinction correction 
is applied. Large filled squares are the results for s tack ed spectra, 
color-coded as for individual galaxies (see section 14.11 1. Big blue 
circles indicate the moment at which models present NUV-r = 4.0 
colors, with diagonal straight lines representing equidistance from 
consecutive models. Also shown on the top right is the median 
error for each individual galaxy. 

showing the fraction of galaxies in each 7 bin as defined 
by the dividing lines in Figure |9] Error bars represent 
fluctuations in each bin, when randomly shifting indi- 
vidual D„(4000) and Hs,a measurements according to 
their corresponding errors; the error bar in each 7 bin 
is then the standard deviation in number counts after 
1000 Monte Carlo simulations. Since 7 correlates with 
quenching speed, that means galaxies with higher 7 val- 
ues will spend less time in the green valley, and are less 
likely to be observed. We show the fractions corrected for 
this (weighted by dy/dt) as dashed lines. We represent 
the data in this work as circles, and the values for low 
rcdshift as triangles. The main conclusion we draw from 
this exercise is that quenching timescales are shorter at 
higher redshift, since the amount of galaxies with higher 
7 is increased. In quantitative terms, this represents a 
factor 2 — 3 decrease in typical transitioning timescales. 

We determined average D„(4000) and Us, a values for 
each magnitude bin by averaging over all galaxies in a 
given bin, weighted by the error in each case. This yields 
an average 7 in each bin, which in turn corresponds to 
a period of time required to cross over the color range 
covered by the green valley. We combine median galaxy 
masses for ea ch bin, as dete r mined from the K-band mea- 
surements by iBundv et al.l (|2007t ) , and number densities 
to calculate a mass flux for each given magnitude. We 
show these values in Table[3] We repeat the procedure in 
the case of the extinction-corrected green valley galaxies, 
and show the results in Table |4j 

Finally, the total mass flux density is the sum of the 



Fig. 10. — Fraction of green valley galaxies as a function of 7 
bins, as shown in Figure(9] All values shown here are for extinction 
corrected samples; circles represent our sample, while the triangles 
are the values found at redshift z ~ 0.1 in Paper I. The blue dashed 
lines in both cases are weighted by dy/dt to correct for the fact that 
galaxies that are quenching faster are less likely to be observed. We 
notice an evolution in the fraction of high-7 values, in that at higher 
redshift the timescales for star formation quenching were shorter. 

mass flux through all magnitudes. We show the re- 
sults in Figure 1111 The upper circle represents the 
value before extinction correction (logp = —0.36 ± 0.08 
M Q yr _1 Mpc~ 3 ), while the bottom value represents 
the mass flux while accounting for extinction correction 
(log/j ex t = -0.80 ±0.51 M yr" 1 Mpc" 3 ). The horizon- 
tal error bars represent the redshift range covered by our 
sample, while vertical error bars include errors in number 
density and transition timescales as described in previous 
sections. 

The first thing to notice is the difference (about a fac- 
tor of 3) between both measurements. This can be at- 
tributed to a combination of factors: on one hand, mea- 
sured number densities are smaller when correction for 
extinction is applied (section 13. 4|) ; on the other hand, 
the contamination of star-forming galaxies also biases 
our sample to higher < Z?„(4000) >, < H$.a > values, 
which in turn represent higher 7's and shorter evolution 
timescales. 

We compare our values to those found in Paper I at 
z = 0.1; the top triangle in Figure [11] represents the flux 
without taking into account dust extinction, and the bot- 
tom one is the extinction-corrected measurement. The 
evolution with redshift from z = 0.1 to 0.8 is evident, 
with mass flux values at intermediate redshifts being 3 
to 5 times higher than those found in the low-redshift 
universe. This reflects the significant change occurring 
in galaxy evolution over cosmic time. 

In Figure [TT] we also show the density growth of the 
red sequence, in units of Mq yr _1 Mpc -3 . This has been 
calculated as follows: the B-band luminosity density at 
the local u niverse has been determined as Jb = 10 7 ' 7 
L© Mpc" 3 (|Bell et al.1120031 iMadgwick et al.l l2002). The 
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TABLE 3 

Mass flux results 
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6.95 
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10.9 
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1.30 


4.18 
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-21.25 
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10.7 


15 


1.32 


4.43 
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10.5 
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TABLE 4 

Mass flux results (corrected for extinction) 
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mass-to-light ratio of red sequence galaxies in the B 
band, in turn, is estimated at (M/Lb) = 6 M /L , 
using a constant, average value for the redshif t range 
covered (jGebhardt et al.ll2003HFaber et alJl2007h . Com- 
bined with the evolution in the number de nsity of galax- 
ies in the red sequence ()Faber et al.ll2007l) . we can then 
calculate the mass flux at each redshift interval as 

/ M \ /l0 A1 °g*o -1\ 

p=3B Kr B ){ iqa^o j^). (*) 

where A log $o is the logarithmic change in the normal- 
ization of the luminosity function at each redshift, and 
Ai(Az) is the change in the age of the universe at each 
redshift interval. 

We find that the mass growth rate was indeed higher 
by a factor of ~ 5 in the earlier universe, pointing to 
another manifestation of downsizing in galaxy evolution. 
Furthermore, our results agree within errors with the es- 
timated growth of the red sequence at z ~ 0.8, reinforc- 
ing the idea that the mass flux through the green valley 
occurred more rapidly in the past. 

4.1. Stacked spectra 

In order to optimize the S/N of the spectral features 
our results are based on, we have produced coadded 
spectra within multiple absolute magnitude bins and re- 
peated the same analysis we applied to our individual 
spectra. The magnitude bins, each lmag in width, were 
defined to be representative of the sample: M_ 24j _ 23 , 
M_23,-22, M_22,-2i, M_2i,-20)M_2o,-i9- We combine 



all individual spectra in our sample within each abso- 
lute magnitude bin performing a meanclip averaging for 
each wavelength pixel between 3800 and 4200 A. The 
resulting stacked spectra are shown in Figure 1121 with 
the shaded region indicating the standard deviation in 
each pixel. With a markedly higher S/N, the spectral 
features in our coadded spectra are clearly distinguish- 
able. To ensure comparison on equal footing with results 
based on individual spectra, we also produced stacked 
specra for absolute magnitudes after applying a cor- 
rection for extinction: M_ 2 5,-24, M_ 2 4,-23, M_ 2 3,-22, 

M_22,-21jM_21,-20- 

Based on the stacked spectra we measure mass flux 
within each magnitude bin repeating the same procedure 
applied to the individual spectra. The measured values 
of D„(4000) and H^ ^ are shown in Figure [9] as filled 
squares. With the improvement in S/N, the results based 
on stacked spectra show a substantial increase in agree- 
ment with the models, strengthening our earlier results. 
We confirm that green valley galaxies show lower values 
for D„(4000) and stronger absorption in H^, indicating a 
stellar population younger than that of their low-redshift 
counterparts, suggesting faster migration from the blue 
cloud to the red clump at higher redshifts. 

The sum of the mass flux through the green valley 
based on the stacked spectra analysis before and after 
extinction is shown in Figure QT] as empty circles, follow- 
ing the same color code as before. The resulting values 
agree with those measured for the average of individ- 
ual galaxies within 1-ct, confirming our results: the mass 
flux density across the green valley at redshift z ~ 0.8 is 
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Fig. 11. — Mass flux density in the green valley. The circles 
represent our data; solid error bars indicate the value calculated 
without extinction correction, while the red symbols with dashed 
error bars include extinction correction. Empty ci rcles indicate 
values measured from stacked spectra (see section 14. It . Values 
at z = 0.1 calculated in Paper I are shown as triangles, and are 
determined with (red) and without (black) extinction correction. 
The solid line sh ows the density gro wth rate of the red sequence, as 
determined from Fabcr ct al. (2007). Our data points show a clear 
increase in the mass flux density across the green valley towards 
earlier times, in agreement with estimates from the growth of the 
red sequence. 

greater than found at low-redshift (z ~ 0.1) by a factor 

Of r-j 5. 

The sum of the mass flux through the green valley 
measured this way is shown in Figure [TT] as hollow cir- 
cles, color-coded as before. The values agree with those 
measured for the average of individual galaxies within 
1-er, confirming our results: the mass flux density across 
the green valley at redshift z ~ 0.8 is greater than found 
at low-redshift (z ~ 0.1) by a factor of ~ 5. 

5. DISCUSSION 
5.1. The build-up of the red sequence 

In iFaber et"al~l ((2007), the authors discuss how the 
red sequence at z — is assembled, from the brightest, 
most massive objects down to the smallest and faintest. 
From DEEP2 and COMB017 data out to z ~ 1, the 
authors are able to measure the evolution of the lumi- 
nosity function and luminosity density - for all galaxies 
and for the blue and red sequences individually - to dis- 
cuss how galaxies might evolve along the color-magnitude 
diagram. The authors argue that the red sequence is 
formed through a combination of star formation quench- 
ing in star forming galaxies and dry mergers, which move 
galaxies along the red sequence towards the brighter end. 

How are the brightest red galaxies formed? One might 
think that the simple quenching of individual blue galax- 
ies, be it by mergers, AGN activity or any other mech- 
anism, would be sufficient to produce the red sequence. 
However, the most massive blue galaxies in the local uni- 
verse are not as massive as the most massive red ones. 




4000 

A. (A) 



4100 4200 



Fig. 12. — Stacked spectra for each magnitude bin. Grey shades 
indicate the standard deviation in each pixel. In all cases spec- 
tral features are clearly distinguishable, in particular the measured 
absorption in at 4100 A. 

Even major mergers, which would produce a galaxy up 
to twice as massive as each individual object, are not 
enough to explain the observed mass function. In fact, 
in the morphological study by Mendez et al. 2011 of 
^300 galaxies in the optically-defined green valley at 
0.4 < z < 1.2, they find that most green galaxies cannot 
be classified as mergers and that, the merger fraction in 
the green valley is in fact lower than in the blue sequence. 

Our data offers an interesting insight into the problem 
of red galaxy formation. Comparison with the mass-flux 
at low-redshift (Paper I) shows that the mass flux oc- 
curred at brighter magnitudes at high redshift (Figure 
[T3|) . indicating the build-up of the most massive end of 
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the red sequence at earlier times, which is in qualitative 
agreement with the evolution of the luminosity function. 
Our sample reliably covers magnitudes down to r ^ — 20 
(not accounting for extinction correction) and r ~ —22 
(accounting for extinction correction) , allowing us to ver- 
ify the decrease in number density and mass flux towards 
fainter objects (see Figure [5]) • This represents a downsiz- 
ing of the green valley evolution, with the red sequence 
forming "from the top down": in the past, more mas- 
sive star-forming galaxies were being formed and sub- 
sequently quenched, forming the more massive red se- 
quence galaxies. At later times, star formation shifts to 
less massive galaxies; these are then quenched as well, 
and the fainter end of the red sequence is created. 

It has been argued before that the massive galaxies in 
the red sequence f orm in earlier times. As an example, 
iBorch et al.1 (|2006f ) have measured the stellar mass func- 
tion of the red sequence up to redshift z ~ 1.0 and have 
noticed that the number density of less massive galax- 
ies has grown more than the number density of massive 
galaxies since that epoch. The evolution we observe in 
the luminosity function supports this idea, with a change 
in the faind-end slope that indicates a deficit of faint, less 
massive galaxies when compared to the same population 
at lower redshift. The analysis of objects in the green 
valley, on the other hand, describes this scenario for the 
first time based on observations of the mass flux of inter- 
mediate galaxies alone. 

Ho wever, this is not the complete picture. (jFaber et al.1 
2007) call upon dry mergers to explain observed proper- 
ties of the brightest red galaxies. These mergers corre- 
spond to the interaction between two or more red se- 
quence galaxies, with little gas involved. The absence 
of copious gas diminishes the subsequent burst in star 
formation that is otherwise expected in wet mergers, be- 
tween gas-rich galaxies. These dry mergers hence result 
in an increase in stellar mass, with no coupled burst in 
activity. It would be interesting to compare the mass flux 
density at the very massive end of the color-magnitude 
diagram to infer whether that is sufficient to create all 
massive galaxies or if dry mergers are a necessity, but 
this is currently an impossible task due to large uncer- 
tainties and small number statistics (especially at these 
limiting magnitudes). 

It could be argued that changes in galaxy properties 
related to an increase in the mass-to-light ratio could ac- 
count for the dimming of red sequence galaxies, in such 
a way that galaxies of a given stellar mass would ap- 
pear fainter at recent epochs. While it is true that the 
mass-to-light ra tio increases as a function of time (e.g. 
iBell et all 120041 ). the change in the r-band is not large 
enough to account for the evolution of the luminosity 
function by itself, especially for passive galaxies. We have 
measured the mass-to- light ratio of galaxies of a constant 
color 4.5 < NUV — r < 5.0 as a function of redshift, 
to evaluate any possible evolution within the color range 
covered by the lumin osity function s, base d on the K-band 
measur ements from iBundv etafl (|2006t l. The average 
value between 0.1 < z < 0.3 is (M/L r ) = 1.9 ±0.6, while 
at higher redshifts, 0.8 < z < 1.0, (M/L r ) = 1.8 ± 0.6. 
We can conclude, then, that the increase in mass-to-light 
ratios of red galaxies is not high enough to account for 
the evolution in the luminosity function, and galaxies in 
that color range are indeed more massive at higher red- 



shift. 

The mechanism through which quenching occurs is 
as yet unclear. As we have seen above, merger activ- 
ity may not be as relevant; AGN activity is also called 
upon, and it is indeed found that green valley galax- 
ies show an increase in AGN fraction - although ac- 
tual correlation with star formation quenching has not 
been unequivocally sho wn (Paper I; iNandra et al.l 120071 : 
ISchawinski et al.l [2009) . In a future paper, we plan to 
investigate this correlation by comparing our results for 
quenching timescales with different tracers of AGN ac- 
tivity, such as X-ray luminosities and optical emission 
line ratios. 
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Fig. 13. — Mass flux density as a function of r magnitude. The 
flux at intermediate redshifts is shown as solid symbols; red tri- 
angles for the extinction-corrected sample, black squares, no ex- 
tinction correction. The lines represents the values found for the 
local universe in Paper I; the red dashed line indicates the extinc- 
tion corrected sample at z = 0.1. In both cases, the peak at high 
redshift is shifted towards brighter magnitudes by > 1 mag. 



5.2. The star formation history of green valley galaxies 

It is important to consider, first, that most of the stel- 
lar mass in red sequence galaxi es is not formed w hile the 
galaxy is in the red sequence (ISalim et al.ll2005h . This 
lends support to our assumption (1), which states that 
galaxies are only moving redward in the color-magnitude 
diagram. However, this is not strictly true; we do know 
a fraction of elliptical galaxies, even in the local uni- 
verse , show signs of recent starburst events (|Thilker et al.l 
2010). If we take into account that a fraction of galax- 
ies in the green valley are turning blue, instead, then the 
inferred mass flux would be smaller. 

At the same time, in Paper I we have argued how 
our choice of star formation histories influences our re- 
sults. In general, assuming shorter periods of constant 
star formation or different models for the time evolution 
decreases the final calculated value for p. Therefore, we 
conservatively assume our measured mass flux density as 
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an upper value for the actual flux in transition galaxies 
from the blue to the red sequence. Still, a decrease by 
0.5 dex brings us to closer agreement with the inferred 
growth of the red sequence. At the same time, the differ- 
ence in final results combined with the observed change 
in the luminosity function and the distribution of 7 val- 
ues makes it safe to assume the evolution in mass flux 
density from z ~ 0.8 to z ~ 0.1 is real. 

Finally, one can argue the star formation history as- 
sumed in this work is simplistic. It has been even argued 
that star forming galaxies follow an inverted r model, 
with an exponentially increasing star formation history 
(|Maraston et al.ll2010D . In any case, a more realistic ap- 
proach to modeling the star formation histories of galax- 
ies could yield different results (likely smaller fluxes, as 
argued above). Thus, we are currently in the process of 
producing models that do not rely on ad hoc star forma- 
tion histories, but instead are drawn from cosmological 
N-body simulations and semi-analytic models, generat- 
ing more physically motivated star formation histories. 
We will compare our spectra with these more refined 
models in a future paper (Martin et al., in prep.). 

6. SUMMARY 

We have shown results of a spectroscopic survey of over 
100 galaxies in the green valley, i.e. transition objects 
with intermediate colors between the blue and the red se- 
quences that are believed to be currently quenching their 
star formation. The data represent the deepest spectra of 
green valley galaxies ever obtained at intermediate red- 
shifts [z ~ 0.8). 

By using measurements of spectral indices - namely 
the break at 4000 A and the equivalent width of the 
H5 absorption line, we are able to infer the star for- 
mation histories of these objects, follo wing the method 
first presented in iMartin et al.l (|2007l) . If one assumes 
an exponential decline in star formation rates, it is then 
possible to measure timescales involved in transitioning 
from the population of star-forming galaxies into the red 
sequence. Combined with independent measurements 
for the number densities in the green valley and typi- 
cal galaxy masses, we are able to measure the mass flux 
density between 0.55 < z < 0.9. 

Our measurements have shown that the mass flux den- 
sity at these redshifts is higher by a factor of ~5 than 
in the low-redshift universe (z ~ 0.1). This can be at- 
tributed to two factors: first, the evolution of the lumi- 
nosity function over the course of the last 6 Gyr, which 
means at z ~ 0.8 the flux is dominated by the brighter 
and more massive objects. Second, quenching timescales 
are shorter at higher redshift, and galaxies transition 
more rapidly from the blue to the red sequence. In addi- 
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tion, we also show that our results are in good agreement 
with estimates for the build-up of mass in the red se- 
quence since 2 = 1. We argue for a scenario in which the 
red-sequence is built "from the top-down" , meaning that 
the most massive objects were quenched at earlier times, 
shifting towards the evolution of less massive galaxies 
today. 

Since we have calculated the quenching timescales of 
individual objects, the next logical step is to compare 
our results with other observable properties in the green 
valley galaxies, in particular merger signatures and AGN 
activity. We expect this analysis will help clarify the role 
of each of these processes in the total quenching of star 
formation, both at low and at high redshifts. 
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